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Abstract: The behavior of simply pulsed qubits (quantum systems with 
two linearly independent states) may be characterized by the energy differ- 
ence AE between the two states of the qubit and by an external stimulating 
potential V{t) that causes transitions between them. Thus, the operation of 
such quantum mechanical systems may be categorized in various regions that 
explicitly depend on AE and V{t). Limiting cases of degenerate, perturba- 
tive, and adiabatic regions are discussed. A comprehensive and illustrative 
map for simply pulsed qubits is presented that can be used as a visual tool 
for students. Furthermore, analytic solutions may be obtained when the in- 
teraction V{t) is proportional to 6{t — tk), namely when a fast interaction, 
called a kick, is used. 



1 Introduction 



Just as classical information and computation are built upon on-off bits, 
quantum information and computation may be constructed on the basis of 
two-state (on-off) quantum systems known as qubits PQE]- These quantum 
building blocks are used to describe and to control atomic and molecular 
reactions [3 E] , electron dynamics j3] , and time ordering in quantum sys- 
tems [6]. However, properties of dynamic two-state quantum systems, such 
as transition rates from one state to the other, are usually found numerically 
since analjd^ic solutions are relatively rare. In this paper we describe some 
simply pulsed qubits where the occupation probability of each of the two 
states, i.e., the likelihood of being in the on state or the off state, can be 
easily controlled. 

The difference between classical and quantum bits is that the classical 
occupation probability of the on state must be either 1 or (definitely on or 
definitely off) at all times, while for a qubit the probability P{t) of being on 
may be any number between and 1 until an observation is made. That is, 
before an observation the qubit may be in a combination or superposition of 
on and off states, similar to a superposition of classical waves. After a qubit 
is observed, P{t) collapses to 1 or for that qubit. If qubits (where is 
large), initially all in the off state, are switched by the same external potential 
and then observed at time t, N ■ P{t) of them will be found to be in the on 
state and A^ ■ (1 — P{t)) will remain in the off state. It is this feature, namely 
that before an observation a qubit can be simultaneously on and off, that 
distinguishes the qubit from the classical bit that simply switches between 
on or off {never being in a superposition of on and off states). In this regard, 
the mathematical rules (i.e., logic) differ for qubits and bits. The more 
sophisticated qubits can in principle perform more complex operations ^1 [7] . 

2 Qubits 

Since bits can only be on or off (e.g., corresponding to small magnets 
that can be magnetized in one of two directions), they are simpler than 
qubits. Before the measurement, a qubit can be in both of its two states 
at the same time. Reliably manipulating this coherent linear superposition 
of states, e.g., with external fields, is essential for quantum computation [T] 
and quantum control jS], but often difficult to achieve. The phenomenon 
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of quantum parallelism |7 allows a quantum processor to execute actions 
in parallel (performing many classical computations at the same time). The 
classical computer, no matter how fast, can only execute commands in series, 
which limits its computational power. Hence, a quantum computer is more 
suitable for solving certain computational problems such as factorization of 
large numbers and database search. Indeed, computer algorithms specifically 
designed to take advantage of the phenomenon of quantum parallelism in 
qubits have been developed |Hlini- Various physical systems that are good 
candidates for qubits are hsted at the end of Section 2. 

2.1 Conceptual Description 

Here we describe a conceptually simple way of controlling a qubit, that 
retains the key features of manipulating more general qubits. The simplicity 
of the following description helps to clarify the concepts without losing the 
essence of more complicated quantum systems. 

We consider a two-state system with possible energies Ei and E2, where 
huj = AE = E2 — El is the energy splitting between these two states. Since 
the zero point of overall energy is arbitrary, the two eigenergies may be 
written as +AE/2 and —AE/2 without loss of generality. The unperturbed 
Hamiltonian Hq does not couple the two eigenstates, i.e., Hq does not cause 
transition of population between the two states of the qubit. Such a system by 
itself is not very interesting since the observable properties do not evolve with 
time. For example, if the system is initially populated with any population 
number A^^i = ■ Pi in state 1 and population A^2 = ^ ■ -f2 in state 2, then 
without an external switching potential A^^i and N2 do not change with time. 

What makes a qubit interesting and useful is the ability to change the 
populations of its states in a controlled way. This is desirable, for example, 
in a device such as a computer, where one wishes to switch the state of a 
storage location back and forth from on to off at will. In classical computers, 
the logic of these binary operations is simple, and is described in terms of 
Boolean algebra ^D], which is a formal way of working with states that may 
be described as either on or off. In such a classical device, the on or off state 
is usually implemented physically by an on or off magnetic field at a given 
storage location. 

Quantum information is similar, except that the system is generally in a 
linear combination (linear superposition) of on and off states. In both classi- 
cal bits and qubits, the switching between on and off states is done with an 
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external field V. Since the state is to be switched from on to off at various 
times, V is some function of time t, i.e., V = V{t). It is this external switch- 
ing field V{t) that enables a computer to manipulate information, performing 
logic operations dynamically. In our model we choose a V{t) that simply cou- 
ples the two states without changing the eigenenergies +AE/2 and —AE/2 
of the unperturbed Hamiltonian Hq. While there are more complex ways to 
manipulate a two-level system, this is sufficient to illustrate how a simple 
qubit works. 



2.2 Mathematical Description 



In quantum mechanics, the state of a system at any time t is fully de- 
scribed by its wave function ip{t). If there are two linearly independent states 



of the system, labeled by 



1 




(e.g., corresponding to on) and 




1 



(e.g., 



corresponding to off), then in the most general case ip{t) is a linear superpo- 
sition of these two states. Mathematically this corresponds to 



Here ai{t) and 02 (t) are the probability (wave- like) amplitudes such that the 
probability of finding ip{t) in the on state is Pi(t) = \ai{t)\'^, and the probabil- 
ity that ilj{t) is in the off state is P2{t) = |a2(^)P- This is analogous to linear 
superpositions of electromagnetic fields in classical wave mechanics, where 
the observed intensity is proportional to the square of either the electric or 
the magnetic field. Normalization condition in quantum mechanics requires 
the constraint Pi{t) + P2{t) — 1 (conservation of probability), meaning that 
no dissipation^ has occurred since the system was initially formed. 

A two-state system can be coupled (i.e., the particle population can be 
transferred between the two states) by applying an external potential V{t). 
In this paper, the external potential has the shape of a single pulse that 
can be sensibly characterized by a single time duration, r. A simple pulse 
could be for example a Gaussian, an instantaneous kick (delta function), 

"'^Dissipation corresponds to leakage of some population out of the two- state system. 
While dissipation can be a significant problem in some practical applications, in this paper 
we assume it is negligible. 



+ a2{t) 




1 



ai{t) 
a2{t) 



(1) 
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or a rectangular pulse. The energy of the system is described by the full 
Hamiltonian H{t), which includes the energy difference term AE and a time- 
dependent external potential term V{t). Here we take AE to be constant in 
time, corresponding to most applications. Hence all the time dependence in 
the system's Hamiltonian H{t) comes from V{t). This full Hamiltonian can 
be written in terms of the Pauli spin matrices (problem 1), namely 

H{t) = Ho + V{t) 

' -AE/2 
+AE/2 

= —cr, + V{t)a,, 

where the widely used Pauli spin matrices [TT] are defined as 
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V{t) 
V{t) 



(2) 
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(3) 



Inserting the Hamiltonian H{t) from expression Q and the wave function 
^p{t) from expression into the time-dependent Schrodinger equation, 



(4) 



we obtain (problem 2) a differential equation for the on and off probability 
amplitudes, ai(t) and 02 (t), as functions of the energy splitting and the 
external potential: 



ih 



di{t) 




_ 02 (t) _ 





-AE/2 V{t) 
V{t) +AE/2 



ai{t) 
02 (t) 



(5) 



Separating the equations, we can also write. 

.^daAt) . , , AE 
iTi — - — = in ai[t) = 



ih 



dt 
da2(t) 
dt 



ai{t) + V{t)a2{t) 
AE 

ih d2{t) = V{t)ai{t)-\ a2(t) 



(6) 
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Conceptually and mathematically, it is convenient to isolate the time 
dependence of the system by defining the evolution operator jT2j or matrix 
U{t,to) that connects the initial state ip{to) with the final state ip{t), 

m = U{t,to)4j{to) . (7) 

All the time dependence of the system is contained in the time evolution 
operator (or Green's function) ?7(t, to), while the initial conditions are spec- 
ified in ip{to). Throughout the rest of this paper, we take to = and write 
the evolution operator as f/(t,to) = U{t). Then, the evolution operator in 
equation ((7j), has the following matrix representation 



f/(t) 



Unit, to) Uu{t,to) 



Unit) U,2{t) 

U2l{t) U22{t) 



(8) 



It is then straightforward to express the time-dependent probability ampli- 
tudes, ai(t) and a2(t), using the evolution operator U(t), namely 



ai(t) 
02 (t) 



U{t) 



ai(0) 
a2(0) 



Unit) Ui2{t) 

U2l(t) U22{t) 



ai(0) 

02(0) 



(9) 



The evolution operator f/(t) may be obtained by solving a differential equa- 
tion very similar to (0): 



Unit) Ui2{t) 

U2l{t) f/22(t) 



-AE/2 V{t) 
V{t) +AE/2 



Unit) t/i2(t) 

f/2l(t) U22{t) 



(10) 



and we require the time evolution operator to start out as the identity matrix 

' Un{0) UuiO) 

f/2l(0) f/22(0) 





" 1 


" 


= I = 





1 


to = 0. 






j;iven by 







(11) 



Unit) f/i2(t) 

f/2l(t) f/22(t) 
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where T is the Dyson time ordering symbol fSJE] that enforces H{ti) acting 
on the system prior to -^(^2) if ti < h- 

Equations jH), (0), ©, ©, ©, dlDl), and (HH) form the basic equations 
for a qubit. These equations are governed by the energy sphtting AE and the 
external potential V{t). Depending on the complexity of the time evolution 
operator, these equations may or may not have analytic solutions. In the 
latter case, numerical solutions of the coupled differential equations (jHI) may 
obscure information about these quantum systems. Analytic solutions, when 
possible, are more convenient and easy to analyze. 

2.3 Practical Considerations for a System of Qubits 

Up to this point, we have treated the qubit as an abstract mathemat- 
ical structure. However, in applications some practical issues have to be 
addressed. For a system of qubits to be capable of performing quantum 
computation [HinilZI, control of population transfer between the two states 
of the qubit is not sufficient. The following list, dubbed the "DiVincenzo 
checklist" [1^] after its developer in 1997, addresses some practical issues for 
a quantum computer: 

• A large number of coupled qubits need to be reliably controlled. 

• It should be possible to prepare the qubits in a desired initial state. 
Commonly, in the initial state all qubits are on or all qubits are off. 

• Decoherence should be minimized. In terms of the energy splitting AE 
between the two states, the decoherence time td over which quantum 
phase information is lost must satisfy td ^ 27rh/AE 

• Quantum gates need to be designed that will control the operation of 
the qubits (these gates play a role similar to classical gates). 

• The information contained in the qubits must be extractable at the end 
of the computation if the outcome of the qubit operations is to be used 
in a productive manner. 

A number of quantum systems satisfy these requirements. Optical pho- 
tons jl6^, nuclear spins [I?! UHl CH] , ion traps [20], nuclear magnetic resonance 
(NMR) [21], and electrons on superfiuid helium [221 some commonly dis- 
cussed examples. 
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The QIST Quantum Computation Roadmap [7] addresses the key issues 
in quantum computation and the latest advances in the field. In the last 
few years, considerable progress has been made in overcoming the problem 
of decoherence, i.e., loss of phase control. Decoherence is an important but 
complicated issue and is discussed elsewhere |2S]- 

Next, we discuss different physical limits for simply pulsed qubits and the 
graphical qubit map is introduced. 

3 Map 

3.1 Significance of the Qubit Map 

A qubit map, such as that shown in Figure 1, is a tool that enables 
one to visualize how the behavior of simply pulsed qubits may depend on the 
variables AE and V{t). In order to make the qubit map a helpful visual tool, 
the map coordinates are taken to be AET/2h (where r is the time duration 
of V{t)) and /q°° V{t')dt' /h. These two variables are dimensionless and range 
in absolute value from to oo. The variable AET/2h determines the effect 
of the unperturbed Hamiltonian Hq on the evolution operator U{t), while 
Io° V{t')dt' /h measures the influence of the external potential on U{t). It is 
useful to think of these variables as independent phase angles or action-like 
integrals ^Hll^. 

It should be noted here that the boundaries separating the qubit map 
regions are not precisely deflned. The regions overlap each other (see curly 
lines in the map) when either of the two phases AEt/2% or /q°° V{t')dt'/h is of 
the order of 27t. This behavior should not be surprising since the deflnitions 
of the degenerate, perturbative, and adiabatic limits are of an approximate 
nature. The map helps to identify the regions where qubit behavior may 
be described analytically, and indicate ballpark values of the parameters for 
which various analytic solutions are applicable. For instance, if we wish to 
investigate slow weakly perturbed qubits, we look for large phases /S.Et/2Ti 
and small V{t')dt' /h. For strongly perturbed degenerate qubits, we look 
for small AET/2h but large values of /q°^ V{t')dt' /%. 

The arrow pointing toward the lower left corner of the Figure indicates the 
"fast" or short-time limit r — 0. Here the time r of the pulse is too short to 
allow for signiflcant phase accumulation due either to the external potential 
V{t) or to the energy splitting AE, i.e., both AET/2h and /o~ V{t')dt'/h are 
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small. The opposite r ^ oo or "slow" limit indicated by the arrow pointing 
toward the upper right corner is where the pulse is long enough for both the 
energy splitting and the external potential to have a large effect over the 
duration of the pulse. 



3.2 Regions of the Map 

It is sensible to separate our map into four regions including the interest- 
ing central region and to examine each one individually. Analytic solutions^ 
exist only for some regions of the map. It should be emphasized that this map 
is vahd only for singly pulsed qubits. Multiply pulsed qubits [211, including 
periodically pulsed qubits j2H] , are not discussed here since an additional time 
parameter, e.g., the typical time between pulses or the period of oscillation 
of the potential, would require a more complicated map. 



3.2.1 Degenerate Region 

When the energy splitting between the two states of the qubits is much 
less than the strength of the external interaction (i.e., AE ^ V{t)), we are 
in the degenerate region of the map. The diagonal extending from the lower 
left ( "fast" ) corner to the upper right ( "slow" ) corner of the Figure forms the 
boundary of the degenerate region: the behavior of the system is dominated 
by the external interaction V{t) whenever we are to the right of or below this 
diagonal line. The evolution of the system is then given by simple sine and 
cosine functions, and the mathematical complexity is greatly reduced. It can 
easily be shown (problems 3 and 4) that as AE — >• in equation (fT^ . the 
evolution matrix of the degenerate qubit or dit^ becomes 



cos a —ism a 
—i sin a cos a 



(13) 



where 



a 



IoVit')dt' 

h 



(14) 



^While solutions of two-state systems with rapidly changing external potentials were 
already examined about 40 years ago the focus has now shifted to quantum control. 
■^Our usage of the term "dit" is not universal. 
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Weakly Perturbative Fields Strong Non-Perturbative Fields 




27r oo 

h 

Fi gure 1: Qubit map for qubits interacting with simply pulsed finite external potentials. 
Here r is the duration of the pulse, A£' is the energy difference between the two states of 
the qubit, and V{t) is the external potential. Note that the single pulse V{t) is neither 
harmonic nor periodic. The two axes are dimensionless. These axes may be thought of 
as action-like |2H1 12^1 phase accumulations due respectively to the energy difference AE' 
between the two possible qubit states and to the external potential V(t). Everywhere 
outside the central region, approximate analytic solutions may be obtained. When the 
total phase associated with the external potential is small, i.e. V{t')dt' /h <C 2tt, then 
the expression (|12|l for the time evolution operator may be expanded in powers of V using 
g-iV'(i )dt /?i ^ _ iv(^t')dt' /h and only the first few terms retained. This corresponds 
to standard perturbation theory for either quickly perturbed qubits where r is small or 
slowly perturbed qubits where r is large. Similarly, if the phase associated with the 
energy splitting is small, AET/2h <C 27r, then we can treat the qubit as degenerate, and 
the solution can be expanded in powers of AE. If both phases are large, then the adiabatic 
approximation generally applies. 
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is the action integral. When a is an integer multiple of 7r/2, the degenerate 
bit becomes a classical bit where the occupation probabilities simply get 
switched between 1 and 0. 

The analytic expression ()13|) is only valid when two conditions I\ETj2Ti <C 
jQV{t')dt'/h and {AET/2hy < jQV{t')dt'/h are simultaneously satisfied. 
The need for the second condition may be removed, thereby extending the 
solution to the entire degenerate region to the right of the diagonal, by re- 
placing a with 

o^iim^, (15) 

n 



where fi(t') = ^V^{t') + {AE/2f. We notice that n{t') = V{t') and thus 
^ = a for ^ 0. 

We will see below that the degenerate region overlaps other regions of the 
map, such as the perturbative, adiabatic, and kicked regions. In these areas 
of overlap, two or more analytic solutions become mathematically equivalent. 



3.2.2 Perturbative Region 



When the external potential V{t) is small, we are in the perturbative 
region. The perturbative region includes both upper and lower left quadrants 
of the qubit map. In this particular region /q°° V{t')dt'/h <^ 27t. The external 
interaction V{t) is sufficiently weak that the qubit remains essentially in its 
initial on or off state, acquiring only an overall phase associated with the 
energy ±AE/2. It can be shown (problem 5a) that the evolution operator 
in equation (fT^ becomes 



f^pert(^) 



-i(t-2t')AE , 

2n' V{t')dt'/h 



i{t-2t')AE 



v{t')dt'/h 

it{AE/2h) 



(16) 



For the extreme special case where V{t) = at all times(problem 5b), the 
evolution operator in equation becomes 



U{t) ^ t/,ero(t) 



^it{AE/2n) g 



(17) 



Setting V(t) = uncouples the differential equations in expression (0) making 
the math trivial. We obtain ai{t) = ai{0)e'^'^^/'^^ and a2{t) = a2{0)e-'^'^^/'^^, 
i.e. the amplitudes of being in the on or off state acquire phases as time 
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evolves but the probabilities |ai(t)p and |a2(^)P remain what they were at 
the initial time t = 0. 

The perturbative region is of little interest in our case since the system 
does not change appreciably from its original state. In the case of zero 
external potential, the qubit does not change its state at all. 

3.2.3 Adiabatic Region 

The adiabatic region covers the upper right, upper left, and lower right 
quadrants of the qubit map, excluding only part of the lower left ("fast") 
quadrant and the central region. The explicit condition for adiabaticity is |21] 

nV{t)AE < [V\t) + (AE/2)2]3/2 , (18) 

which can be further simplified in the degenerate and nondegenerate limits. 
In the degenerate adiabatic case AE <C V{t), so (fTS|) yields the condition 
hV(t)AE < V^{t) or HAE/t < V^it). The time derivative of the external 
potential is of the order l^(t)/r, i.e., V ~ l^/r, which is the typical slope of the 
single-pulse potential. In the nondegenerate adiabatic case AE ^ V{t), and 
the validity condition ^ becomes hV{t) -C {AEf or hV{t)/T < {AEf. 

The meaning of the adiabatic region is that the external potential changes 
slowly, so that the quantum system is able continuously to adjust to the new 
Hamiltonian. Consider for example a classical system of two pendulums 
connected by a spring. The system can oscillate in one of two normal modes: 
a lower-frequency swinging mode where the two pendulums move left and 
then right in unison, and a higher-frequency vibrational mode where the 
spring is alternately stretched and compressed. If we move the pivot point of 
one or both pendulums slowly and carefully enough, the pendulums continue 
to oscillate in the same (lower or higher) mode they started in, because not 
enough energy is provided by the external perturbation to switch the motion 
from the lower-frequency oscillation to the higher-frequency oscillation or vice 
versa. Similarly, an adiabatically driven qubit that is initialized in the lower 
energy {—AE/2 or off) state before the pulse is turned on will at any future 
time remain in the lower energy state of the full Hamiltonian H{t). After the 
pulse is over and the Hamiltonian returns to the unperturbed Hamiltonian 
Hq, the qubit must return to the off state with energy —AE/2. Thus, the 
adiabatic region is of no great interest in the case of pulsed qubits since we 
wish to completely transfer the population of a qubit from one state to the 
other. 
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Mathematically, the evolution operator for an adiabatically driven qubit 

is 



COS 9{t)+^^^ sin 9{t) -^^sm9{t) 

-i^^ sin 9{t) cos 9{t) sin 9{t) \ ' ^ ' 

where 9{t) and Q{t) are defined as in Section 13.2. H and for simplicity we 
assume V{t) = V{0), i.e., the external potential returns to its initial value at 
the end of the pulse. Equation (fT^ can be reduced (problem 6) to equations 
we presented above for the degenerate and zero external potential cases. 



^adiab(^) — 



3.2.4 Central Region 

None of the limits discussed previously overlap with the central region 
of the map. Here the two independent variables AET/2h and V{t')dt' /ft 
are both of order 27r, and there is no large or small parameter to simplify 
the solution of the problem. If the pulse has a particularly simple (e.g., 
rectangular) shape, an analytic expression is possible even in this region, 
but for a general pulse shape no analytic solution exists and the evolution 
must be computed numerically. 



4 Kicked Qubit Approximation 

There is a useful approximation for simply pulsed qubits, namely the 
sharp, narrow pulse or "kick" limit in which the width r of the pulse goes to 
zero, while the integrated strength or the area under the external potential 
curve 

ak= [\{t')dt'/h (20) 

remains fixed. Formally, the shape of a very narrow pulse of finite total 
strength a may be expressed by a delta function: V{t') = akh6{t — tk), where 
tk is the time at which the pulse is applied. The kicked region corresponds 
to the lower half of the qubit map in Figure 1. Here the duration of the 
pulse is so short that AET/2h <^ 27r, i.e., there is not enough time for the 
energy splitting AE to have a significant effect while the pulse is active. 
The integrated strength of the pulse, a^, may be either large or small in 
this region. If is large, we are in the lower right quadrant of the map. 
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where the kicked region overlaps with the adiabatic region. If is small, 
we are in the lower left quadrant, where the kicked region overlaps with the 
perturbative region. 

In the kicked region, the evolution of the qubit can be described as follows: 
(i) before the kick, the qubit evolves in accordance with the unperturbed 
Hamiltonian Hq, i.e. the on and off probability amplitudes ai{t) and 02 (t) 
acquire phases proportional to the energy splitting AE (as in (jl7|) ). while the 
on and off probabilities, given by |ai(t)p and |a2(t)P, remain unchanged, (ii) 
For the duration of the kick, the energy splitting associated with the unper- 
turbed Hamiltonian Hq may be ignored, and a transfer of population between 
on and off states may occur, depending only on the integrated strength ak 
of the kick, (iii) Finally, free evolution governed by Hq resumes after the 
kick is complete, and no further population transfer occurs. This sequence of 
events is very similar to the collision approximation studied in introductory 
physics, where (i) a particle initially moves freely with constant velocity, (ii) 
then undergoes an instantaneous collision with a wall, during which the ve- 
locity is changed but there is not enough time for the particle to move, and 
(iii) finally, the particle once again resumes free flight with a (new) constant 
velocity. 

For such kicks, the integration over time in equation ()12|) is straightfor- 
ward (problem 7) and the time evolution matrix becomes a product of three 
parts. 



t^kicked(t) = e*2ft^* ^'''"^e •''k-' e'2n^k^- 

^iAEt/2h _-^iAE(t-2t,)/2h ^-^ 

_^^-iAEit-2t,)/2n ^-^AEt/2n ^^g 



(21) 



Here, we may set^ the Dyson time ordering symbol T 1 for the duration 
of the pulse [t^ — e < t < + e). It can also be shown (problem 8) that 
equation (PT|) will reduce to other limiting cases discussed previously. 
Moreover, the occupation probabilities for a kicked qubit initially in state 1 
are, from equation Q, 



Pi{t) = \ai{t)\^ = \Un kicked (i^)P = cos^ftfc 
P2{t) = |a2(t )|' = \U2i kicked (t) I' = sin^ak • (22) 

^For a single kick setting T ^ 1 corresponds to the limit of no time ordering, which 
is beyond the scope of this paper. For some useful discussions on time ordering see 
reference [HI- 



14 



Since sin^ + cos^ = 1 for this closed two-state system, the conservation 
of population holds, i.e., there is no dissipation. We notice also that the 
transfer probability P2{t) does not depend on the energy splitting AE. This 
simple instructional example may be extended to a series of kicks [21] • It is 
one of the few cases in which analytic solutions may be obtained for qubits 
controlled by external potentials. 

5 Summary 

A primary motivation of this paper has been to give an overview of simply 
pulsed qubits. We have developed a simple, yet comprehensive and instruc- 
tive, two-dimensional map for such qubits. Analytic solutions with limited 
applicability are found and their corresponding regions on the map are iden- 
tified. One promising class of analytic solutions for dynamic qubits are the 
kicked solutions, valid for qubits subject to fast interactions. These pro- 
vide clear and useful examples to students and offer an alternative to the 
well-estabhshed rotating wave approximation (RWA) j2Hll2niEni- The RWA 
method assumes a sinusoidal external perturbation with frequency chosen to 
be in resonance with the desired transition between two states. The solutions 
presented in this paper are, in principle, valid for a wide spectrum of pulse 
widths T and energy splittings AE. 

Qubits are building blocks for quantum computation and quantum infor- 
mation. Their behavior and interaction with each other and the environment 
have to be addressed before a quantum computer becomes a reahty. Ulti- 
mately, N-qubit systems and their interactions will have to be successfully 
controlled. The simply pulsed qubits discussed in this paper may form a 
basic building block for more complex interconnected N-qubit systems. 
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APPENDIX 
Problems for Students 



Problem 1: (a) Show that a 2x2 matrix that corresponds to a two- 
state quantum mechanical observable can be written in terms of the 
Pauli matrices. [Hint: Write down a generic Hermitian 2x2 matrix 
and try to expand it as a sum of the Pauli matrices with appropriate 

prefactors]. (b) If the qubit is implemented as a spin-1/2 particle, ^ 


1 

represents a state where the spin points in the —z direction, then what 
do operators a^, ay, and correspond to physically? 



represents a state where the spin points in the +z direction, and 



Problem 2: Using equations ((T]), (j2)), and (jl)), verify equation 
Rules for simple matrix operations are available in many elementary 
math books and on the web. 

Problem 3: (a) Solve the evolution of the degenerate qubit or dit 
by plugging AE = into equation (0), thus proving that the prob- 
ability amplitudes, ai(t) and a2(t), oscillate sinusoidally with a = 
jQV{t')dt' /h. (b) Make plots of the probabilities Pi = |ai(t)p and 
P2 = |a2(^)P- What can you deduce from the plots about dissipation 
in the system? (c) Find the condition for complete transfer from state 
1 to state 2. 

Problem 4: (a) Solve the evolution of the degenerate qubit or dit by 
using the evolution operator in equation (fT^ and then plugging U{t) 
into equation Q. [Hint: Since AE = 0, Hq vanishes and time ordering 

effects disappear (T = 1). Then U(t) = e~^Jo ^W'^^/^^ ^^d the following 

r Q 

identities may be useful: e*^'^^ = I cos{6) + zcr^ sin(^^) = ^ 

and e*^°"^ = /cos(6') + ia^ sin(^^).] (b) Do you get the same probability 
amplitudes as in problem (3b)? (c) Check to see that the evolution 
operator is unitary, that is U{tyU{t) = I. 

Problem 5: (a) Show that in the perturbative limit, the evolution 
operator in equation (fTB|) satisfies equation (fTUj). (b) Do the same for 
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the case of zero external potential. That is, prove equation (fT7|) . You 
may use the result from part (a). 

• Problem 6: (a) Show that the adiabatic equation (|19|) for the evolution 
operator reduces to the degenerate equation (fT^ in the limit /S.E 0. 
(b) Show that the adiabatic equation (|19p reduces to the zero potential 
equation (fT7|) in the limit V(t) 0. 

• Problem 7: Prove equation (j2H) by integrating and then multiplying 
out the exponentials in the first line, thus obtaining the matrix in the 
second line of equation (j?T|) [Hint: You may find useful the identities 
given in problem (4)]. 

• Problem 8: (a) Show that equation (PT|) for a kicked qubit is equiva- 
lent to the perturbative equation (fTBj) when the kick strength is small, 
i.e., afc <^ 1 [Hint: When <C 1, cosa^ ~ 1 and sina^ ~ a^. In equa- 
tion replace t' with since we are in the kicked region in this 
problem], (b) Show that the equation for a kicked qubit reduces 
to the degenerate equation (fT^ in the limit AE 0. 
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